MASSACHUSETTS INSTITUTE OF TECHNOLOGY 
ARTIFICIAL INTELLIGENCE LABORATORY 



A.L Memo 498 October, 1978 



CALCULATING THE REFLECTANCE MAP 



/""N 



Berthold K. P. Horn 
Robert W. Sjoberg 



Abstract 

It appears that the development of machine vision may benefit from a detailed 
understanding of the imaging process* The reflectance map, showing scene radiance as a 
function of surface gradient, has proved to be helpful in this endeavor. The reflectance 
map depends both on the nature of the surface layers of the objects being imaged and the 
distribution of light sources. Recently, a unified approach to the specification of surface 
reflectance in terms of both incident and reflected beam geometry has been proposed. 
The reflecting properties of a surface are specified in terms of the bidirectional 
reflectance-distribution function (BRDF). 

Here we derive the reflectance map in terms of the BRDF and the distribution of 
source radiance. A number of special cases of practical importance are developed in 
detail. The significance of this approach to the understanding of image formation is 
briefly indicated. 
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,-s §1. THE REFLECTANCE MAP 

The apparent "brightness" of a surface patch depends on the orientation of the 
patch relative to the viewer and the light sources. Different surface elements of a non- 
planar object will reflect different amounts of light towards an observer as a 
consequence of their differing attitude in space. A smooth opaque object will thus give 
rise to a shaded image, one in which brightness varies spatially, even though the object 
may be illuminated evenly and covered by a uniform surface layer. This shading 
provides important information about the object's shape and has been exploited in 
machine vision [1-8]. 

A convenient representation for the relevant information is the "reflectance map" [4, 
6]. The reflectance map, R(p , q), gives scene radiance as a function of surface gradient 
(p f q) in a viewer-centered coordinate system. If z is the elevation of the surface above 
a reference plane lying perpendicular to the optical axis of the imaging system, and if x 
and y are distances in this plane measured parallel to orthogonal coordinate axes in the 
image, then p and q are the first partial derivatives of z with respect to x and y. 



f*\. . 



p = dz/dx and q = dz/dy 

The reflectance map is usually depicted as a series of contours of constant scene 
radiance (Fig. 1)- It can be measured directly using a goniometer-mounted sample, or 
indirectly from the image of an object of known shape. Alternatively, a reflectance 
map may be calculated if properties of the surface material and the distribution of light 
sources are given. One purpose of this paper is to provide a systematic approach to this 
latter endeavor. Another is to derive the relationship between scene radiance and image 
irradiance in an imaging system. This is relevant to machine vision since gray-levels are 
quantized measurements of image irradiance. 



§2. MICROSTRUCTURE OF SURFACES 

When a ray of light strikes the surface of an object it may be absorbed, transmitted, 
or reflected. If the surface is flat and the underlying material homogeneous, the 
reflected ray will lie in the plane formed by the incident ray and the surface normal 
and will make an angle with the local normal equal to the angle between the incident 
ray and the local normal. This is referred to as "specular," "metallic," or "dielectric" 
reflection. Objects with surfaces of this kind form virtual images of surrounding 
objects. 

Many surfaces are not perfectly flat on a microscopic scale and thus scatter parallel 
incident rays into a variety of directions (Fig. 2a). If deviations of the local surface 
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FIGURE 1: A typical reflectance map for a surface, with both a glossy and a matte 
component of reflection, illuminated by a point source. The coordinates are surface 
slope in the x and y directions, and the curves shown are contours of constant scene 
radiance. 
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normals from the average are small, most of the rays will lie near the direction for ideal 
specular reflection and contribute to a surface "shine" or "gloss" 

Other surface layers are not homogeneous on a microscopic scale, and thus scatter 
light rays which penetrate the surface by refraction and reflection at boundaries 
between regions with differing refractive indices (Fig. 2b). Scattered rays may re- 
emerge near the point of entry with a variety of directions and so contribute to 
"diffuse," "flat," or "matte" reflection. Snow and layers of white paint are examples of 
surfaces with this kind of behavior. Frequently both effects occur in surface layers, 
with some rays reflected at the nearly flat outer surface of the object, while others 
penetrate deeper and re-emerge after multiple refractions and reflections in the inhomo- 
geneous interior. 

The distribution of reflected light in each case above depends on the direction of 
incident rays and the details of the microstructure of the surface layer. Naturally, what 
constitutes microstructure depends on one's point of view- Surface structures not 
resolved in a particular imaging situation are taken here to be microstructure. When 
viewing the moon through a telescope, for example, smaller "hillocks" and "craterlets" 
are part of this microstructure. This consideration leads to more complicated models of 
interaction of light with surfaces than those discussed so far. It is possible, for instance, 
to consider an undulating surface covered with a material which in itself already has 
complicated reflecting behavior (Fig. 2c). 

Reflectance is not altered by rotating a surface patch about its normal when there is 
no asymmetry or preferred direction to either the pattern of surface undulations or the 
distribution of sub-surface inhomogeneities. Many surface layers behave this way and 
permit a certain degree of simplification of the analysis. Exceptions are such things as 
diffraction gratings, iridescent plumage and the mineral called "tiger eye " These all 
have a distinct directionality in their surface microstructure and will not be considered 
here any further. 

Considerable attention has been paid to the reflective properties of various surface 
layers. Some researchers have concentrated on the experimental determination of 
surface reflectance properties [9-21]. At the same time, many models have been 
developed for surface layers based on some of the considerations presented above 
[22-35]. Models often are too simple to be realistic, or too complicated to yield 
solutions in closed form. In the latter case, Monte-Carlo methods can be helpful, 
although they lead only to numerical specification of the reflecting behavior. Purely 
phenomenoiogical models of reflectance have found favor in the computer graphics 
community [36, 37, 38]. Several books have appeared describing the uses of reflectance 
measurements in determining basic optical properties of the materials involved [39, 40, 
41]. Attention has been paid, too, to the problem of making precise the definitions of 
reflectance and related concepts [42, 43]. 
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§3. RADIOMETRY 

A modern, precise nomenclature for radiometric terms has been promoted by a 
recent NBS publication [43]. The following short table gives the terms, preferred 
symbols and unit dimensions of the radiometric concepts we will have occasion to use 
for the development presented here: 



RADIOMETRIC CONCEPTS 

Radiant flux # [W] 

Radiant Intensity / = d$/dia [W'sr -1 ] 

Irradiance E = d$/dA [W • m -2 ] 

Radiant Exitance M = d$/dA [W»nr 2 ] 

Radiance L = d 2 $ / (dA • cos $ • dw) [W.m- 2 »sH] 

Radiant flux, $, is the power propagated as optical electromagnetic radiation and is 
measured in watts [W]. The radiant intensity, /, of a source is the exitant flux per unit 
solid angle and is measured in watts per steradian [W»sH]. The total flux emitted by 
a source is the integral of radiant intensity over the full sphere of possible directions 
(47r steradians). The irradiance, E, is the incident flux density, while radiant exitance, 
Af, is the exitant flux density, both measured in watts per square meter of surface 
[W • m -2 ]. The total radiant exitance equals the total irradiance if the surface reflects 
all incident light, absorbing and transmitting none. 

The radiance, L, is the flux emitted per unit foreshortened surface area per unit solid 
angle. Radiance is measured in watts per square meter per steradian [W • m -2 • sr -1 ]. It 
can equivalently be defined as the flux emitted per unit surface area per unit projected 
solid angle. Radiance is an important concept since the apparent brightness of a surface 
patch is related to its radiance. Specifically, image irradiance will be shown to be 
proportional to scene radiance. 

Radiance is a directional quantity. If the angle between the surface normal and the 
direction of exitant radiation is 0, then the term "foreshortened area" stands for the 
actual surface area times the cosine of this angle 6. Similarly the "projected solid angle" 
stands for the actual solid angle times the cosine of the angle 6. Here we will use the 
symbol w to denote a solid angle, while SI will be used to denote a projected solid angle. 
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If du> and dQ are corresponding infinitesimal solid angles and projected solid angles 
respectively, then ^ 

dQ = dw * cos 8 

The following example (Fig. 3) will illustrate some of these ideas. Consider a source 
of radiation with intensity / in the direction of a surface patch of area dA, oriented 
with its surface normal making angle $ with the line connecting the patch to the source. 
In fact, as seen from the source, it appears only as large as a patch of area dA-cos 
would when oriented perpendicular to this line. The corresponding solid angle is simply 
the area of this equivalent patch divided by the square of the distance from the source 
to the patch. Thus, 



r\ 



du> = dA * cos 6 / r 2 [sr] 
The flux intercepted then is 

</$=/. do = / . dA • cos 6 / r 2 [W] 

The irradiance of the surface is just the incident flux divided by the area of the surface 
patch. 

E = d$/dA = / . cos 6 / A [W. m~2] 



(T\ 



SURFACE 
NORMAL 




FIGURE 3: Point source illuminating a surface, 
illustrating basic radiometric concepts. 
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§4* THE BIDIRECTIONAL REFLECTANCE-DISTRIBUTION FUNCTION 

The Bidirectional Reflectance-Distribution Function (BRDF) was recently introduced 
by Nicodemus, Richmond, Hsia, Ginsberg, and Limperis [43] as a unified notation for 
the specification of reflectance in terms of both incident- and reflected-beam geometry. 
The BRDF is denoted by the symbol f r and captures the information about how 
"bright" a surface will appear, viewed from a given direction, when it is illuminated 
from another given direction. To be more precise, it is the ratio of reflected radiance 
dL r in the direction towards the viewer to the irradiance dEf in the direction towards 
a portion of the source. In symbols, 

f r Wi><l>i*r><t>r) = ^M >*&'*?&) J <»/(*/. 0/) [sH] 

Here, and together indicate a direction, the subscript / denoting quantities 
associated with incident radiant flux, while the subscript r indicates quantities associated 
with reflected radiant flux [43]. 

The geometry is as depicted in the figure (Fig. 4). A surface-specific coordinate 
system is erected with one axis along the local normal to the surface and another 
defining an arbitrary reference direction in the local tangent plane. Directions are 
specified by polar angle 6 (colatitude) measured from the local normal and azimuth 
fs angle (longitude) measured anti-clockwise from the reference direction in the surface. 

In general, incident flux may arrive from many portions of extended sources, so incident 
radiance £/(0,-,0/) is a function of direction. If we consider the component of flux d$f 
arriving on the surface patch of area dA from an infinitesimal solid angle dwj in the 
direction (tf/ f 0/) we obtain 

</#/ = Lj • cos 6j • du>i • dA = L x ♦ dQ t • dA = dE t • dA 

where dEf—L^oosd^dooj is the incident irradiance contributed by the portion of the 
source found in the solid angle dw { in the direction (0/,0/). Similarly, the radiant flux 
emitted into an infinitesimal solid angle dw r in the direction (0 r ,0 r ): 

(ft$ r = dL r • cos 6 r • du> r • dA = dL r • dQ r • dA 

where dL r {8 rf <t> r ) is the radiance in the direction (# r ,0 r ) due to the reflection of the 
flux incident from direction (0/ f 0/). The notation dX % where X is one of the 
radiometric quantities introduced in §3, will always denote a directional quantity, that is, 
one which depends on either the incident or exitant direction. The notation d 2 X will 
mean a bidirectional quantity, which depends on both the incident and exitant 
directions. Thus, the incident flux </$,• depends only on the direction of incidence, but 
^^ the exitant flux d 2 ^ r depends on both the direction of emittance and (implicitly) on the 

direction of incidence. 
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FIGURE 4: Local geometry of incident and 
reflected rays needed for the definition of the 
bidirectional reflectance-distribution function 
(BRDF). Redrawn from [43]. 



From these values of incident and exitant flux, the BRDF is defined as follows: 

/r(0/.0/: O r <<t>r) - i^r/Mr) / <»/ = dL r /dE t 

and thus has dimension inverse steradian [sH]. The BRDF allows one to obtain 
reflectance for any defined incident and reflected ray geometry simply by integrating 
over the specified solid angles [43]. 



§5. INTEGRALS OVER SOLID ANGLES AND PROJECTED SOLID ANGLES 

The admitting aperture of an imaging system may occupy a significant solid angle 
when seen from the point of view of the objects being imaged. We will furthermore 
have to deal with extended sources. In both cases it is necessary to integrate various 
quantities over solid angles or projected solid angles. This can be accomplished by 
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double-integration with respect to the polar and azimuth angles (Fig. 5). If AT is the 
quantity to be integrated, we have 

J^X dw = J^Jq X sin 6 dd d<t> 



and 



Jr r« r*/2 

^ X dQ. = J ^ J o X cos sin d$ d<f> 



If for example X — 1 and the region of integration is the hemisphere above the object's 
surface, then 

w Z^ = J_J sin d6 d* - 2ir 



while 



% fr /*tt/2 



J^ X </fi = jTj[* (1/2) sin 20 dd d<f> = n 
The latter result will be used in the discussion of perfectly diffuse reflectance. 



FIGURE 5 : Polar and azimuth angles used in 
double integrals over specified solid angles. 
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§6. PERFECTLY DIFFUSE REFLECTANCE 

A perfectly diffuse or "lambertian" surface appears equally bright from all directions, 
regardless of how it is irradiated, and reflects all incident light [43]. Thus the reflected 
radiance is isotropic, that is, L r is constant, with the same value for all directions 
(6 r ,<t> r ). Also, the integral of reflected radiance over the hemisphere above the surface 
must equal the irradiance E. This implies that the BRDF for this ideal surface, f ridt is 
constant, and that the radiant exitance M equals the irradiance E. If the reflected 
radiance is L r , then the radiant exitance can be found by integration: 

M aa L L r dQ. r — L r ir 

From this one finds that 

f rM - L r /B, - L r /M = l/w 



f~\ 



If we have an extended source with radiance £,-, then the irradiance on the surface due 
to a small portion of solid angle de>; lying in the direction (0/,0/) is 

dEi(fi h 0,) = L,($ it ty) cos $i dw { . 

The reflected radiance is then 

L r xm (1/7T) j^ LjiBi, 0,) COS 0, ; dW; 

This is a form of Lambert's cosine law. 

§7. COLLIMATED SOURCES AND THE DIRAC DELTA-FUNCTION 

Not all sources are extended. One way to deal with sources that are highly 
collimated is to treat them as limiting cases of extended sources, with the distribution 
tending towards an impulse or delta-function. If this is to be expressed in a coordinate 
system of polar and azimuth angles, one has to take into account the non-uniform 
spacing of coordinates. Consider a collimated source which produces an irradiance E 
on a surface oriented orthogonally to the direction (0 O , O ) of its rays. Clearly the 
radiance L t of this source should be zero except for this direction. The product of 
Dirac delta functions, 8(6,- - O ).5(0 ; .- O ), will be a useful ingredient of the 
formula expressing L t as a function of the angles. One must insure, however, that the 
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Horn & Sjoberg 13 October, 1978 

irradiance on a surface lying orthogonal to the rays equals E& that is, 

nw/2 
L t sin $ t d6i d<t>f. 

Clearly this can be accomplished if 

L t = E 8($i - 6 ) 6(0/ -' 06) / «n 6 

This is called the "double-delta" representation of source radiance for a coliimated 
source. It can also be written in an alternate form using the identity 

6 [/(*) " /M - * <* - *tf / ■fW- 
where f'(x ) is the derivative of f(x) evaluated at x = Xq, provided that f(x) has an 
inverse in the domain of interest. (This identity can be confirmed by integrating each 
side with respect to x over the relevant interval.) Then, 

L t = E «(cos $i — cos 6 ) 8(0/ — O ). 



§8* PERFECTLY SPECULAR REFLECTANCE 

A perfectly specular or "mirror-like" surface reflects light rays in such a way that 
the exitant angle 6 r equals the incident angle 0/ and that the incident and reflected ray 
lie in a plane containing the surface normal. The reflected radiance of a surface patch 
in the direction (0 r ,0 r ) is simply the source radiance in the corresponding reflected 
direction. That is, 

L r (6 r ,<f) r ) = L/(0 r ,0 r + n) 

The surface thus forms a virtual image of the source. From the definition of the 
BRDF, we see that 



L r = J/ r dE, - J^ f T L t dSl t 



That is, 

% tt pr/2 



r^ 



nv/2 
f r L; cos 9 t sin $i d6i dfy 

We can satisfy the conditions stated above if we let 

fr.is = *( e i ~ Q r) *(0/ - 0r + *) / ( sin .*/ ^ AO 
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This is called the "double-delta" form of the BRDF for perfectly specular reflectance. 
Using the identity mentioned in the last section, we can write this in an alternate form 
[43]: 

fr.is =26 (sin20 r - sin20,.) b (0 r - 0,- + n) 

§9. ANALYSIS OF IMAGE-FORMING SYSTEM 

We will now analyze a simple image-forming system (Fig. 6). We assume that the 
device is properly focused; that is, those rays originating from a particular point on the 
object which pass through the entrance aperture are deflected to meet at a single point 
in the image plane. Similarly, rays originating in the infinitesimal area dA$ on the 
object's surface are projected into some area dA p in the image plane and no rays from 
other portions of the object's surface will reach this area of the image. Further, we 
assume that there is no "vignetting," that is, the entrance aperture is a constant circle 
of diameter d and does not become occluded for directions which make a larger angle 
with the optical axis. The effect of vignetting on image irradiance will be considered 
later. 

The exposure of film in a camera is proportional '^\ai^^i&ii^'E p " % 2s^'^t^- 
levels in a digital imaging system are quantized measurements of image irradiance. In 
order to calculate image irradiance we must first determine the flux d$ L passing 
through the entrance aperture arriving from the patch of area dAq on the object. 

d *L = d4)Jsi L r d$l r 

where Q. r is the projected solid angle subtended by the aperture. We will also need to 
know the area dA p of the image of the patch, since image irradiance E p is the flux per 
unit area: ■ 

E p = d$ L /dA p 

If 0/ is the angle between the normal on the surface and the line to the entrance 
aperture nodal point, while a is the angle between this line and the optical axis, then, 
by equating solid angles, 



Consequently, 



f 0mm S i 



(d^ cos 0/)// o 2 _ (dA p cos a)// 2 



E p = (fo/fp) 1 cos a j^ L r (cos r /cos $/) du> r 
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FIGURE 6 : A simple image forming system. Light collected by the lens from 
the surface patch of area dA§ is projected into the image patch of area dA p . 



^**N 



Here the integral is over the solid angle occupied by the entrance aperture as seen from 
the patch on the surface. If we assume that the lens is small relative to its distance 
from the object, 6 r is approximately the same as 0/, and the ratio of their cosines is 
unity. Furthermore, the reflected radiance L r will then tend to be constant and can be 
removed from inside the integral. The solid angle occupied by the lens as seen from the 
surface patch is approximately equal to the foreshortened area (n/4) d 2 cos a, divided 
by the distance (/q/cos a) squared. Finally, one obtains the well-known result, 

E p = L r (n/4)(d/f p )2cos4a 

That is, image irradiance is proportional to scene radiance. The factor of 
proportionality is n divided by four times the effective f-number (f p /d) squared, times 
the fourth power of the cosine of the off-axis angle, a. Thus the "sensitivity" of such an 
imaging system is not uniform over an image, but is constant for a particular point in 
the image. Vignetting introduces an additional variation with image position. Ideally, 
an imaging device should be calibrated so that this variation in sensitivity as a function 
of a can be removed. 
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Other kinds of imaging systems, such as microscopes or mechanical scanners, lead to 
somewhat different expressions. Generally, however, image irradiance is proportional to 
scene radiance in such systems too. At this point we should remember that scene 
radiance depends on properties of the surface layer (BRDF) and the distribution of 
light-sources (source radiance) since 



' r ~ Jfl,- fr L i 



d$l: 



§10. VIEWER-ORIENTED COORDINATE SYSTEM. 

So far we have considered directions from the object to the image-forming system 
and to light sources in terms of a local coordinate system with one axis lined up with 
the surface normal. Such coordinate systems will vary in orientation from place to 
place and are thus inconvenient for the specification of global distributions such as that 
of source radiance. A coordinate system fixed in space will be more suitable, 
particulary if one of the axes is lined up with the optical axis (Fig. 7). In this viewer- 
oriented coordinate system we introduce polar angle 6 measured from the r-axis and 
azimuth angle measured from the *-axis in the plane perpendicular to the *-axis. 
Here, the z-axis is parallel to the optical axis. Directions to sources of light can be 
given using these two angles. If the sources are far away in comparison to the size of 
the objects being imaged, then source radiance will be a fixed function of these angles 
independent of the point on the surface being considered. 




FIGURE 7: Viewer-oriented, global coordinate 
system useful for specification of the distribution 
of source radiance Lj. 
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§11. THE SURFACE NORMAL 

In the local coordinate system the surface normal lies along one of the axes, or 
equivalently, it is the direction corresponding to zero polar angle. In the viewer-oriented 
coordinate system the surface normal will correspond to some direction, say ($„,<£„). 
The corresponding unit vector is, 



n = (cos 0„ sin 6„ , sin 0„ sin $„ , cos $„ ) 



The surface of the object may be specified by giving "elevation" z as a function of the 
coordinates x and y. We can give an expression for the surface normal in terms of the 
first partial derivatives of z with respect to jc and y, if these exist. Let the first partial 
derivatives be called p and q. Then the vectors (1,0, />) and (0, 1,?) are tangent to the 
surface, as can be seen by considering infinitesimal steps in the x and y direction. The 
surface normal is perpendicular to all vectors in the tangent plane and so is parallel to 
the cross-product of these two: 

(l,0,p)X(0,l,q) = (-p,-g,l) 

Thus the unit normal can be written 

n = (-p,-q,l) / yj\ + p2 + q 2 

The following results are obtained by equating terms in the two expressions for the 
surface normal: 

sin $„ = V/> 2 + ? 2 / V* + P 2 + <P- 
cos $ n = 1/ Vl + P 2 4- q 2 

sin 0„ = -q / V/> 2 + q 2 
cos 0„ = -p / yjpl + ql 



Conversely, 



p — —cos 0„ tan d n 
q = —sin 0„ tan 6„ 



^"\ 
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^ §12. RELATIONSHIP BETWEEN LOCAL AND VIEWER-ORIENTED 

COORDINATE SYSTEMS 

In order to calculate the scene radiance, we will integrate the product of the BRDF 
and the source radiance over all incident directions. Since the BRDF is specified in 
terms of the local coordinate system, while the distribution of source radiance is likely 
to be given in the viewer-oriented coordinate system, it will be necessary to convert 
between the two. Given the direction of the surface normal, (0 n ,0„), and the direction 
to a portion of the source (0,,0,), both specified in the viewer-oriented system (Fig. 8), 
we have to find the incident direction (0/,0,) and the exitant direction (0 r , r ) both 
specified in the local system. Alternatively, given the surface normal and the incident 
direction we may have to find the direction to the source and the exitant direction. 
Note that r = $ n since the exitant ray lies along the z-axis in the direction towards 
the viewer. Further, since we have excluded anisotropic surfaces, we are only interested 
in the difference between r and 0,-. From the relevant spherical triangle (Fig. 9) we 
obtain 



/■"s 



r^ 



A 




FIGURE 8 : Surface normal and direction to portion of the source shown in 
viewer-oriented coordinate system. 
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8 r = 6 n 




/^*\ 



r\. 



FIGURE 9: Spherical triangle extracted from previous figure 
and used in derivation of transformation equations between the 
surface-normal, local coordinate system, and the viewer-oriented, 
global coordinate system. 



Cosine formula: ~ " ~~ 

cos 0j =cos 6 S cos 6 n + sin 6 S sin 6 n cos (0 5 — <f>„) 
Sine formula: 

sin 6 i sin (0 r — 0/) = sin 6 S sin (0 5 — 0„) 
Analogue formula: 

sin e t cos (0 r — 0/) = cos 6 $ sin 6 n — sin 6 S cos n cos (0^ — n ) 

The Jacobian of the transformation from ($ s ,<t> $ ) to (0/,0/) equals 

Mi/des) O0//a0 5 ) - (30//30 5 ) (30,/30 5 ) - (sin 5 /sin tf f ). 

(The Jacobian will be required below when converting a double integral with respect to 
one set of coordinates to one in terms of the other.) The above formulae allow us to 
find the incident direction from the source direction. Quite symmetrically, we can also 
obtain the source direction from the incident direction: 



Cosine formula: 

cos 6 S = cos $j cos 6 r + sin B f sin B r cos (0 r — 0/) 
Sine formula: 

sin 6 S sin (0 5 — n ) = 
Analogue formula: 



sin $j sin (0 r — 0/) 
logue formula: 
sin 6 S cos (0, — 0„) = 



cos ; - 



?/ sin r — sin 0/ cos r cos (0 r — 0/) 



Horn & Sjoberg 20 October, 1978 



fs 



^*\ 



The Jacobian of the transformation from (0 j9 0,) to (6 s ,<t> s ) equals 

Ws/Mi) 0*0, /30/) - (d4> s /d$i) (de s /d<f>i) = (sin 0,/sin 6 S ). 

§13, SCENE RADIANCE 

It follows from the definition of the BRDF that reflected radiance can be written as 
the integral 

Using polar and azimuthal angles this becomes 

nir/2 

Here we integrate over all possible incident directions (0,-,0/) and calculate source 
direction (0 f ,0 5 ) from the given surface normal (0„,<t>„) and the incident direction. 
The inner integral has the limits to ir/2 for 0,-, corresponding to directions within the 
hemisphere visible from the surface. The integration can be extended to the full sphere 
of directions if the integrand is forced to be zero when 0/ lies between ir/2 and ir. This 
can be accomplished by replacing cos 0/ by max [0, cos 0,]. Hence 

L r = J _^J o f r ^ max[0, cos 0,] sin 0/ d6 t dty 

Since the integral now is over the full sphere of directions, it can be easily rewritten 
using any other set of polar and azimuth angles. Using the viewer-oriented coordinate 
system, for example, we obtain 

L r = J v J o f r L t max[0, cos 0,] sin 0, d$ s d<f> s 
That is, 

fM#i'> O r <<t>r) I/(0,,0,)-max[O, cos 0,] sin $ s d6 s d<f> s 

Here we integrate over all possible source directions (6 s ,<t> s ) and calculate incident 
directions (0 z -,0 ; ) from the given surface normal (6„,<j) n ) and the source direction. We 
now have two convenient forms for the calculation of scene radiance. We proceed to 
calculate reflectance maps for a few simple combinations of BRDF and distributions of 
source radiance. 



rs 
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§14a. COLLIMATED SOURCE, LAMBERTIAN REFLECTANCE 

For a lambertian reflector, f r — l/n. For a collimated source, 

L t = E 8 (6 S - 6 ) 6 (0, - o )/sin 6 

where E is the irradiance measured perpendicular to the beam of light arriving from 
source direction (6q,<P ). Substituting into the second form of the expression for scene 
radiance above, we obtain 

(E /ir) 8(6 S - d ) S(0 5 - O ) max[0, cos $ t ] (sin 6 s /sin O ) d0 s d<t> s 

This is equal to 

L r = (E /n) max[0, cos Q(] 

where 6/ = 0/ when 8 S = O anc * $s — #o- * n ^ s case * 

COS $/ =sb COS r COS Q + s * n #r S * n ^0 °° S (#0 "~ $91 ) 

and 

cos (00 — n ) = cos O cos n — sin O sin rt 

To obtain the reflectance map, scene radiance as a function of surface gradient, we can 
substitute expressions in p and q for these trigonometric expressions. The result is 



Rip ^4) — (E /n) max 



f V+PqP+W) 



where 



Pq = —cos O tan O 
go = —sin O tan 6 

A surface with gradient (p ,q ) is normal to the direction of the incident light rays. 



rs 
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§14b. UNIFORM SOURCE, LAMBERTIAN REFLECTANCE 

A uniform source has constant incident radiance. Let L t = Lq. Again, for a 
tambertian reflector, / r = 1/tt. Substituting into the first form of the expression for 
scene radiance, we obtain 

rn/2 

This becomes 



nir/2 
(Lq/k) cos $i sin dj d$i dty 



J rn/2 
[ sin 20, -^ := £<, 



Not surprisingly, the reflected radiance is independent of the surface orientation in this 
case. 



§14c. HEMISPHERICAL UNIFORM SOURCE, LAMBERTIAN 
REFLECTANCE 

A hemispherical uniform source is described by 

Li(6 s ,<t> s ) = Lq for e s < tt/2 
Li(d s ,<f> s ) = for $ s > it/2 

In order to evaluate the double integral for scene radiance, it is necessary to know the 
value 0/ of the incident angle 0,-, which corresponds to the. "horizon* of the "sky,* that 
is, $ s = ?r/2. From the coordinate transformation equations one can easily see that 

cot 0/ = —tan r cos (0 r — 0,) 

For r — 7r/2 < 0,- < r + jt/2, the horizon cutoff only occurs for $( > ir/2 and 
can be ignored, since the inner integral is from 0, = to 0,- = tr/2 only. For the 
other half of the range of /t this cutoff must be considered. Now, 

nir/Z 
f r L { cos $i sin fydOidti 

so, 

nmin[0/,7r/2] 
cos 6; sine; ddidfy 

Letting = 0,- — r and 0' = 0,- — r + tt, we can split the outer integral, 
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L r - (lo/ir) J [ w/2 J q cos 0, sin*,- d6 ( d<f> + (!<,/*) J^ J Q cos 0, sin 0,- d9 t d<f/ 

Now, 

J cos f sin 0,- rf0, = 1/2 
so the first term is simply Lq/2. Next, note that 

J cos ^ sin0,. d$i d$i = (1 - cos 20/)/4 = sin2 0//2 
where, since cot 0/ = —tan r cos (0 r — 0,), 

sin 2 0/ = 1/[1 + tan 2 6 r cos 2 (0 r - <f>#\ 
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FIGURE 10: Cross-section through uniform hemispherical source and surface 
element, illustrating horizon cutoff and portion of extended source not visible 
from surface. 
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The second term thus becomes ~ 



'tt/2 
jr/2 

which equals 



Tf/2 

(L /2n) J 1/[1 + tan* r cos* ^ d<f>' 



f~\ 



(L /2w) [cos r tan-' (cos r tan 0)]^ = (Lq/2) cos r 
Adding up the two terms we finally get 

L r(9n><t>n) = (V 2 ) + Cos 0„) = Lq cos* (0„/2) 

This is the result found by Brooks [44]. From it the reflectance map can be found 
immediately: 

R(p,q)*=(E /2)(l + l/yjl+ p 2 +q 2 ) 

§15a. COLLIMATED SOURCE, SPECULAR REFLECTANCE 

For specular surfaces, 

f r = s($i - r ) 6(0,- - r + ir) / (sin 0,- cos 0,) 

Using the source radiance from §14a above, and the first form of the expression for 
scene radiance, we obtain 

nr/2 
(Lo/sin O ) 5(6; -$ r ) *(0 / _0 r + 7r ) 5(0, -0 O ) 5(0, -0 O ) d$ t dty 

That is, 

L r = Lq 5(0/ - O ) 5(0/ - O ) / sin O 

where 0/ and 0/ are the values of 0, and <p s corresponding to 0,- = r and 
0/ = r + »r. Using the equations for the coordinate transformations, one finds that 
0/ = 20 r and 0/ = 0„. Thus, 

L r = Lq 5(20 r - O ) 5(0„ - O ) / sin O 
and finally, 

L r (d n ,<t>„) = (Iq/2) 5(0„ - o /2) 5(0„ - 0q) / sin O 
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To express this as a function of p and q we have to remember that 

6[f(x,y)-f(x 0t y )] 5[g(* , y) - g(x , y )] = 6(x-x ) 8 (y - y )/J(x , y ) 

where /(a:,^) is the Jacobian of the transformation from (x,y) to (f,g), (provided this 
transformation has an inverse in the region of interest): 

J(x , y) = (df/Sx) (dg/dy) - (df/dy) (dg/dx) 
The Jacobian of the transformation from (p,q) to (0„,<t>„) is 

J(P*1) = l/l V/> 2 + ? 2 + P 2 + ? 2 ) ] 
Let 

/>l = —cos O tan 6 /2 
q x = —sin O tan o /2. 

Then, noting that sin O = 2 sin (0 o /2) cos (0 o /2), one can write 

sin 6 = 2 V/>j 2 + ?i 2 / (l 4- ?i 2 + ?, 2 ) 
and therefore, 

*(/; , 9) = (E /4)5 (p - /?,) 6 (q - ?,) (l + P{ 2 + ^2)2 

Thus, a surface element with gradient (p { ,q x ) is oriented to specularly reflect the 
collimated source towards the viewer. This gradient can be related to the gradient (pq, 
q ) introduced earlier: 

Pi - PQ ( Vl + Pq 2 + q 2 - l) / <j> 2 + q 2 ) 

<7l = % (Vl + Po 2 + ?0 2 - I) / 0>o 2 + ?0 2 ) 

When the point (p ,q ) is not far from the origin, then (/>j,?,) is approximately 
midway between the origin and (p , q ). 



/■"> 
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rs §15b. UNIFORM SOURCE, SPECULAR REFLECTANCE 

It is easy to see that for a specular surface under a uniform source, the scene 
radiance will be constant and equal to the source radiance. 

L r — Lq 

This is the same result as the one we obtained for the uniform source and lambertian 
reflectance. Thus a diffuse surface appears just as bright as a specular surface if both 
are viewed with uniform illumination. In fact, all surfaces reflecting the same fraction, 
P say, of the total incident light will appear equally bright under this illumination 
condition. 



§15c. HEMISPHERICAL UNIFORM SOURCE, SPECULAR REFLECTANCE 

In this case, 



L r (Q n ,<t>„) = Lq for 0„ < jr/4 
£ r (^,0„) = O for 6„> n/A 



The reflectance map is 



R(p,q) = £0 for p 2 + g 2 < 1 

R(P,q) = for/>2 + ql > 1 

SUMMARY AND CONCLUSIONS 

We have shown that image irradiance is proportional to scene radiance and that 
scene radiance depends on surface orientation. The reflectance map gives scene 
radiance as a function of the gradient. It can be calculated from the bidirectional 
reflectance-distribution function (BRDF) and the distribution of source radiance. 
Several special cases were worked out in detail. Each could have been developed more 
easily by a direct method, but was obtained from the general expression for scene 
radiance to illustrate the technique. The general expression allows one to find the 
reflectance map even if the source radiance distribution or the BRDF is only given 
numerically. 
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